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ABSTRACT 


Authors  usually  neglect  the  effect  of  Coulomb  field  in  the 
space  of  radius  e>^(3kT/2)  surrounding  a  charged  particle,  and 
derive  equations  of  the  Fokker-Planck  type  from  the  Liouville 
equation  of  a  plasma.  effect  of  the  inner  core 

field,  neglected  usually,  is  shown  to  be  present  in  an  equation  of 
the  Fokker-Planck  type  as  an  additional  term  which  is  similar  to 
the  collision  term  in  an  equation  of  the  Boltzmann  type.  It  is  shown 
that  the  order  of  this  additional  term  may  easily  be  larger  than  those 
of  the  friction  and  diffusion  terms,  the  effect  of  the  outer  core  field. 
When  the  effect  of  the  outer  core  field  is  larger  than  that  of  the  inner 
core  field,  the  assumption  of  binary  interaction  is  shown  to  be  not 
feasible.  The  possible  nonlinear  and  non-Markovian  behavior  of  a 
particle  is  neglected. 
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LIST  OF  MAIN  SYMBOLS 

- 

A 

friction 

"b 

diffusion 

®a 

electric  charge  of  particle  a 

Ea 

macroscopic  electric  force  on  particle  a 

Fa»  Fj, 

one-particle  distribution  functions 

^ab»  ^ac 

two-particle  distribution  functions 

^ab 

^ab  "  ^a  ^b 

• 

Ld 

Debye  length  (the  radius  of  the  outer  core  field) 

[see  Eq,  (2.4)j 

ma 

mass  of  particle  a 

Na 

total  number  of  particle  a  in  volume  V 

"a 

Na/V 

■Pa 

momentum  of  particle  a 

■^a 

position  vector  of  particle  a 

R 

the  radius  of  the  inner  core  field  ^see  Eq.  (2.5)j 

'} 

time 

t  ^ 

V 

volume  of  the  container 

S 

-L 
n  S 

• 

‘Aab 

potential  energy  between  two  particles  a  and  b 

iv 

potential  outside  the  outer  core 


^abji 


potential  inside  the  outer  core 


potential  inside  the  inner  core 

the  average  period  of  correlation  in  the  outer  core 
field  ^see  Eq.  {3,l)j 

f  the  average  period  of  correlation  in  the  inner  core 

field  ^see  Eq.(3.2)j 
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I.  INTRODUCTION 


We  investigate  the  process  of  deriving  the  equation  of  the  one- 
particle  distribution  function  of  a  plasma  and  show  that  it  is  convenient 
to  give  the  interaction  term  as  the  sum  of  the  term  of  collision  which 
is  characteristic  of  the  Boltzmann  equation  and  those  of  diffusion  and 
friction  which  are  proper  to  the  Fokker-Planck  equation. 

Usually  authors  (Refs.l,  2  and  3)  derive  equations  of  the  Fokker- 
Planck  type  by  assuming  that  the  probability  of  a  pair  of  particles 
interacting  with  each  other  at  distances  shorter  than  e^/(3kT/2)  is 
negligibly  small.  The  interaction  outside  this  distance  (and  inside  the 
Debye  length  distance)  is  weak  and  of  a  fairly  long  correlation  time. 
Therefore,  it  is  possible  to  expand  the  s-body  function  similarly  to  the 
Mayer  cluster  expansion  and  the  interaction  in  this  outer  core  is  con¬ 
veniently  presented  by  the  terms  of  friction  and  of  diffusion.  Here  we 
pose  the  question  whether  it  is  proper  to  neglect  the  effect  of  the  inner 
core  of  the  field. 

According  to  the  result  in  this  paper,  the  effect  of  the  inner 
core  field  may  easily  be  larger  than  those  of  the  friction  and  diffusion 
terms  caused  by  the  outer  core  field.  The  probability  of  the  inner  core 
correlation  is  small.  But  the  correlation  itself  is  strong.  Thus  the 
total  effect  may  not  be  negligible. 
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With  respect  to  the  inner  core  field,  usual  methods  of  expanding 
correlation  functions  similar  to  Mayer's  cluster  expansion  may  not  be 
readily  applied,  as  the  expansion  may  not  converge.  In  this  paper,  we 
present  the  effect  of  strong  interaction  of  short  correlation  periods  in 
the  inner  core  by  "collision"  and  the  weak  interaction  of  long  corre¬ 
lation  periods  in  the  outer  core  by  "friction  and  diffusion."  As  a 
result,  the  equation  of  a  one-body  function  includes  the  "collision" 
term  and  "friction  and  diffusion"  terms.  Unlike  the  case  of  neutral 
particles,  two  particles  immediately  before  a  "collision"  by  means  of 
the  inner  core  field  are  correlated  to  each  other  due  to  the  outer  core 
field.  We  note,  however,  that  weak  correlations  in  the  outer  core  are 
effective  because  of  long  correlation  periods.  These  weak  correlations 
are  not  effective  when  we  integrate  them  over  short  periods  of  inter¬ 
actions  in  the  inner  core.  In  the  short  period  of  interaction  in  the 
inner  core,  only  the  strong  interaction  developed  in  the  inner  core  is 
effective.  Therefore  we  may  neglect  the  history  (or  memory)  formed 
in  the  outer  core,  when  we  consider  the  interaction  in  the  inner  core. 
This  is  our  assumption  which  makes  the  present  investigation  possible. 

In  the  future,  we  may  consider  the  interaction  between  these  two  kinds 
of  correlations  which  are  to  be  one,  but  now  separated  rather  artificially. 
Related  to  this  point  of  the  problem,  we  note  that  Tchen  (Ref.  4)  recently 
investigated  the  interaction  between  the  correlation  caused  by  the  outer 


core  field,  and  the  effect  of  the  macroscopic  field  caused  by  the  macro¬ 
scopic  nonuniformity  of  a  plasma.  In  this  paper,  these  two  kinds  of 
interactions  are  neglected  as  secondary  effects. 


II.  ASSUMPTION  OF  BINARY  CORRELATIONS 


We  consider  a  mixture  of  particles  of  different  kinds  and  assume 
that  all  the  particles  are  electrically  charged  and  that  their  interactions 
are  caused  by  Coulomb  electrostatic  force.  The  distribution  function  of 
a  single  particle  belonging  to  a  group  will  be  denoted  by  F^,  etc. 

The  binary  correlation  function  for  a  pair  of  particles  will  be  denoted 
by  F^tj,  F^^,  etc.  By  integrating  the  Liouville  equation  with  respect 
to  the  position  and  momentum  coordinates  (q,  "p)  of  all  the  particles, 
except  for  one,  we  obtain 


'JL 

dt 


fftQ 


'i 


(2.1) 


Here 


(2.2) 


and  the  summation  X  is  to  be  extended  to  all  the  particles  except  for 
the  one  under  consideration.  For  the  equation  of  F^j,,  we  obtain 
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(2.3) 


( 


A 

dt 


+ 


A.  +  A  .  ^  <Mito 
dTo  '"b  <?q^  ^5^ 


According  to  the  results  of  Tchen's  analysis  (Ref.  1),  we  separate  the 
domain  of  f"^  -  j  in  three: 


Domain  I: 


>  i-o  *  fX 

'  b 


AvN^et)^ 

VkTo 


(2.4) 


Domain  II  (Outer  Core): 


•-0^  l%-Tbl>Rab‘ 

Domain  III  (Inner  Core): 

«<lb^  l^-^bl 


9‘ 

(SkT/2  ) 


(2.5) 


(2.6) 


Here  3kT/2  is  the  average  energy  of  a  particle.  In  Domain  I,  we  may 
write  for  F^jj  simply 

Pot  -  Po  Pt 

The  interaction  of  particles  a  and  b  is  effective  only  when  the 
distribution  is  not  spatially  uniform.  Of  course,  we  neglect  the  weak 
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correlation  introduced  by  the  spatial  nonuniformity  which  was  considered 
by  Tchen  recently  (Ref.  4).  In  Domain  II,  we  write  for 


^o^b-^  ^ab' 

^abc’‘ ^ab~  (2.8) 


Because  of  Eq.  (2.5),  it  is  permissible  to  do  so,  as  carefully  studied 
by  Tchen  (see  Ref.  1,  Appendix  A).  In  Domain  III,  however,  any  rapid 
convergence  of  an  expansion  similar  to  Eq.  (2.8)  may  not  be  expected. 
Separating  the  domain  of  the  integration,  with  respect  to  q^ ,  we  write 
for  Eq.  (2.1) 


Instead  of  separating  the  domain  of  integration  with  respect  to  qj^ ,  it 
might  be  convenient  to  separate  the  potential  into  three  parts  as 

shown  in  Fig.  1  (see  following  page).  The  singularities  at  the  boundaries 
of  domains  are  simple  and  may  not  cause  any  serious  difficulties  in  our 
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mathematical  treatment.  Concerning  Eq.(2.9)>  we  may  define  the 
macroscopic  electric  field  by 

^  *  -ST  f  “^0 

We  write  for  Eq.  (2.9) 
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Our  further  task  is  to  obtain  in  Domain  II  and  in 

Domain  III.  The  terms  on  the  right-hand  side  of  Eq.  (2.3)  may  be 
neglected  if  we  assume  that  the  nonlinear  correlation  is  negligibly 
small.  After  neglecting  these  terms,  there  remains  no  term  of 
integration  with  respect  to  and  Therefore,  it  is  possible  to 

write  Eq.  (2.3),  separately  in  Domain  II  and  III,  as  follows: 

In  Domain  II, 


By  considering  Eq.  (2.1),  we  substitute 


for 


/  d  ,  %  ,  d 

'■jr*  mT  "SC 


in  the  above  equation.  We  neglect  the  nonlinear  terms  which  appear  in 

the  equation  after  the  substitution,  and  also  ^ 

t'Qo  ^Pa  ^Pb 

which,  according  to  Tchen's  study,  is  not  significant  so  far  as 
Domain  II  is  concerned.  Therefore,  for  the  equation  to  which  ^'ab 


is  subject,  we  obtain 

/  d  ,  'Pa  d  ,  d  , 

,  is^.  IfsL.  Fft-  Fa 

Equation  (2.12)  is  valid  only  in  Domain  II. 


(2.12) 
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In  Domain  III,  neglecting  the  nonlinear  terms  in  Eq.(2.3),  we 


have 

/  d  ,  Pa  d  ,  d 

d^ab  d  d4>ba  d  ^  (2.13) 

-''it*'’ 


III.  CORRELATIONS 


In  Domain  II,  the  average  period  of  correlations  is  long,  although 
the  correlations  are  weak. 

T  ^0 

^  *  {3kT/m)!t  (3«1) 


gives  the  order  of  the  average  period.  In  Domain  111,  the  average  period 
of  correlations  is  short  although  the  correlations  are  strong. 


- 


R 

(3kT/m) 


(3.2) 


We  assume  that  the  ratio  of  two  periods  is  sufficiently  larger  than  unity. 


(3.3) 
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Since  the  changes  of  variables  in  Domain  II  during  the  period  Tjg  are 
assumed  to  be  negligibly  small,  we  obtain 


dPb 


dqij  dpb  ds 


(3.4) 


by  averaging  each  term  of  Eq.(Z.ll)  over  rjjf  .  Herej  dashed 
characters  denote  respectively  the  averages  of  the  original  functions 
over  Tjn  along  the  trajectories  of  concerned  particles: 

^  •  T5"  J  Fq  (t  +s)ds 

0 

(i)  The  last  term  of  the  right-hand  side  of  Eq.(3.4)  may  be 
calculated  by  means  of  Eq.  (2.3).  Here  we  simply  mention  the  elaborate 
method  by  Kirkwood  and  Ross  (Ref.  5).  They  consider  a  phase  space 
transformation  function  Substituting 

t+s)»fK**^(q,p/q',p';S)  F^^(q',p';f  )d'q' dp' 

for  Fj^Ij  in  Eq.(2.3),  we  obtain  the  equation  of  We  do  not  repeat 

their  investigation  in  detail.  However,  we  must  remember  that  the 
necessary  assumptions  by  which  the  last  term  of  Eq.  (3.4)  is  reduced 
to  the  collision  integral  are  not  precisely  given  in  our  case.  The 
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correlation  developed  in  Domain  II  does  exist  before  the  two  particles 
begin  to  interact  in  Domain  III.  As  mentioned  previously,  however,  we 
assume  that  the  correlation  in  Domain  II  is  weak  and  that  its  integral 
over  Tjjj  is  not  significant  as  compared  with  the  integral  of  the  strong 
correlation  newly  developed  in  Domain  III.  This  collision  term  may  be 
calculated  by  assuming  that  the  potential  between  two  particles  is  the 
inner  core  potential  as  given  in  Fig.  1.  The  collision  term  thus  obtained 
is  denoted  by 


'  St  4)//.  ^^0 


(3.5) 


(ii)  By  averaging  each  term  of  Eq.(2.12),  over  ,  we  obtain 


(3.6) 

dS,  "  dtb  d$b  “ 


Here,  of  course,  we  assume  that  in  Domain  II,  the  changes  of  functions 
in  T]a  are  small.  By  integrating  the  equation  along  the  trajectories 
of  two  particles,  we  obtain 


Fab(f) 


SFq  i  d<^ba  ^  dFb 
Spa  ^  <^9b  Spb 


(3.7) 


-JO- 


The  integrand  is  considered  to  be  a  function  of  qa{s),  q^Cs), 
Qa^s)=qg(th  (t~s)pg/m^ 


is  neglected  in  Domain  II,  p  and 


Since  .  i-^b-  i^ 

dqo  dPa  dqt  dPt  . . 

are  considered  invariants.  Substituting  Eq.(3.7),  we  obtain  for  the 


first  term  of  the  right-hand  side  of  Eq.{3.4), 

t-rji 


d^b  (s  ) 

Here  — ^22 —  is  neglected  in  the  non-Markovian  approximation, 

oPa 

Defining  A  and  B  by 


Nb  f  d^ab(f)[  / 

lL 


FdFtjds]  dqti  dpb 


(3.8) 


(3.9) 


the  interaction  in  Domain  II  is  represented  by 


—  •AFa-h-i^  -T  Fa 
OPq  oPaPPa 


(3.10) 
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Summarizing  the  above  results,  the  equation  for  is  given  by 


■h  Sq  Eg- 


/d&.) 

'  dt  'coll. 


(3.11) 


In  order  to  calculate  A  and  B,  Tchen  uses  the  Fourier  integral 
representation  of  and  5 


^5(0^  *  «o  expCiu  •q)Y  (u) 

-CD 

/^OO 

dy  »xp( it ‘q)2otJly} 


(3.12) 


Taking  Maxwellian  distributions  for  F^  and  Fjj  in  the  integrands, 
he  investigated  A  and  B  in  detail. 


IV.  EVALUATION  OF  INTERACTIONS 


According  to  Eq.(3.8),  we  see  that  the  order  of  A  is 


(4.1) 


Considering  Eq.  (3.11)  as  the  basic  equation,  we  may  now  examine  the 
relative  effects  of  the  terms  on  the  right-hand  side.  That  of  B  is 
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(4.2) 


°[^h(T^)  iji 


Here  c  is  the  average  speed  of  particles 


C  =  OkT/m) 


he 


(4.3) 


and  we  may  take  for  the  order  of  operator  d/d^ 


“  (3kTm)he 


(4.4) 


Considering 


.  2  _  /<r 

"  47rne2 


.  d^'BFa 


we  obtain  for  the  orders  of  and  ^  . 

dp  dp  dp 


We  may  estimate  the  order  of 

o[(^)  ,.]=7rR^^n= - IT  "VS — T 

L\  dt  /CO//J  3Xff27r)'«m^Lo’ 


(4.5) 


(4.6) 


(4.7) 


Making  ratios  of  these  values,  we  have 

o[(dF^/dt)coii.]  o[(dF^/df)coii]  3 

_L - ^  »47rr.Lo 

o/  J  o[d^T  Fa/dp^j 

We  interpret  Eq.(4.8)  as  follows; 


By  putting 


(4.8) 


(4.9) 
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we  have 


LD/8  - 


(kT)i 
e  n'^ 


(4.10) 


From  Eqs.{2.5)  and  (4.5),  we  have 

Lp  jx 

R  477-)  e’  n 
It  is  easily  shown  that 

or 

109  (-^J  *  log(67r)-h  2  log(^) 

If  we  take  Lq/§>  I ,  it  is  clear  that  8/R  >1  .  When  Lq/8  =  I  , 
we  obtain  from  Eq.  (4. 10) 

(47r)^ 


or 

log,oO  *  3  log,QT  +  5.04 


(4.11) 


(4.12) 


(4.13) 


(4.13a) 


(4.14) 


(4.14a) 


where  n  is  in  cm"^  and  T  in  °  K.  Equation  (4.14a)  is  plotted  in 
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Fig.  2.  If  is  extremely  large,  we  may  neglect  completely  the 

effect  of  microscopic  encounters  and  the  equation  is  of  the  Vlasov  type. 
When  is  fairly  larger  than  unity,  the  interaction  by  the  outer 

core  may  be  neglected  and  the  equation  is  of  the  Boltzmann  type.  As 
^D/8  decreases,  the  friction  and  diffusion  terms  appear  as  subsidiary 
ones.  Near  Ld/8  =  /  ,  the  effect  of  the  inner  core  (collision)  and  the 

effect  of  the  outer  core  (friction  and  diffusion)  may  be  comparable  to 
each  other.  As  Lq/8  decreases  further,  the  collision  term  may  be 
negligible  and  the  equation  is  of  the  Fokker-Planck  type.  At  the  same 
time,  the  assumption  of  binary  encounters  becomes  unplausible.  In 
other  words,  when  the  equation  of  the  Fokker-Planck  type  is  plausible. 


Fig.  2  n  Is  Given  By  Unit  cm"^  and  T  by  °K  In  The  Domain 
Above  The  Line,  ^-D/8  >  1  • 


encounters  are  to  be  multiple  (Ref.  3).  In  Fig.  3,  these  various  domains 


are  indicated  schematically. 


'09io'^ 


Pig.  3  tog  IQ  T  -  log|Qn  Space  («  Divided  By  Lines,  ^0/Z  sconsiant, 

In  Pour  Domains. 

The  Equation  (e: 

(1)  Of  The  Vlasov  Type  Wftien  ^D/B  is  Extremely  Larger 

(2)  Of  The  Boltsmann  Type,  as  ^0^  Decreases; 

(3)  Of  The  DoUtmann  Type  -f  Fokker-Planck  Type 

(Binary),  When  Approaches  To  Unity; 

(4)  Of  The  Fokker>P!anck  Type  (Multiple),  As 
LO/8  Decreases  Further. 

(n  Domain  (4),  Encounters  Are  Multiple  (More  Than  Binary). 

V.  CONCLUDING  REMARKS 


By  the  order  analysis  above,  we  see  that  the  friction  and  diffusion 
terms  due  to  binary  interactions  have  validity  only  as  corrections  to  the 
collision  term  of  the  Boltzmann  equation.  Within  the  limits  of  binary 
interaction  these  Fokker-Planck  terms  by  themselves  appear  to  have  no 
region  of  validity.  When  the  collision  term  of  the  Boltzmann  type 
becomes  negligible  compared  to  the  long-range  interactions,  we  must 
consider  multiple  (more  than  binary)  interactions. 
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culating  (  ^|j  {3.5},  a  particle  which  contributes 

to  this  tern,  is  to  have  a  velocity  satisfying  the  condition  that  the 
nearest  distance  between  these  colliding  particles  is  shorter  than  R. 
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